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Abstract. - The steady states of two vibrated granular gases separated by an adiabatic piston 
are investigated. The system exhibits a non-equilibrium phase transition with an spontaneous 
symmetry breaking. Even if the gases at both sides of the piston have the same number of 
particles and are mechanically identical, their steady volumes and temperatures can be rather 
different. The transition can be explained by a simple kinetic theory model expressing mechanical 
equilibrium and the energy balance occurring in the system. The model predictions are in good 
agreement with molecular dynamics simulation results. The macroscopic description of the steady 
states is discussed, as well as some physical implications of the symmetry breaking. 



The adiabatic piston is a prototypical construction to 
investigate the relaxation of a system and the relevance 
of fluctuations and successive dynamic time scales as the 
system evolves. It consists of a gas-filled container divided 
into two compartments by a freely moving adiabatic pis- 
ton [1]. When the initial states of the gases in the two 
compartments are different, the relaxation can be rather 
complex even when the whole system is isolated and the 
final state is the equilibrium one PHI]. Of course, the 
behavior of the system is much less understood when it 
can not relax to an equilibrium state due to some external 
constrains. In this case, even the answer to apparently 
simple questions is not known. For instance, is there any 
thermodynamic-like parameter characterizing the macro- 
scopic equilibration between the two subsystems? If there 
is a final steady state, does it depend and in which way 
on the properties of the piston? 

Granular media are a key type of systems to address 
the above and many other issues related with fundamen- 
tal concepts in far from equilibrium systems are granular 
media. Due to the inelasticity of collisions, granular fluids 
are intrinsic non-equilibrium systems since no steady state 
is possible in an isolated system. On the other hand, the 
energy dissipation in collisions offers the possibility of new 
energy balances and, consequently, steady states show- 
ing many peculiarities as compared with time- independent 
states of molecular, elastic systems are observed [7l|8]. 

In the case of granular gases, the dynamics is domi- 
nated by (inelastic) binary collisions. Therefore, it is not 
surprising that the techniques developed in kinetic the- 



ory and statistical mechanics for ordinary molecular flu- 
ids have been extended with success to them [9l[T0]. An 
advantage of granular gases over molecular fluids is that 
their smallest scale is macroscopic and, therefore, many 
of their characteristic features are directly observable in 
experiments. 

The aim of the work being reported here is to study the 
steady state of two vibrated granular gases separated by 
an adiabatic piston. There are some previous studies of 
adiabatic pistons involving inelasticity effects. In some of 
them, the whole system is isolated, and the granular gases 
collide elastically with the piston. The system as a whole 
is cooling, and the two granular gases are considered ini- 
tially in the same macroscopic homogenous state. It has 
been shown that this system exhibits a phase transition 
with an spontaneous symmetry breaking [Ullll] . In other 
studies, infinite baths of elastic gases colliding inelastically 
with the piston have been considered. It is assumed that 
the baths are not affected by the piston, so that a steady 
state is reached after some transient [13j[T4]. Of course, 
all the above studies deal with highly idealized situations 
hard to approximate in real experiments. In particular, 
the homogeneous cooling state is of too short a duration 
to be observed. To overcome this limitation, in the present 
work, energy is continuously supplied to the gases in both 
compartments by vibrating the walls. Actually, the con- 
struction described below tries to mimic the experimental 
setup used in ref. [TS] , although the results reported there 
can not be related to the theoretical analysis here, due to 
the rather different range of densities in both cases: very 
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Fig. 1: Sketch of the system considered in this work. All the 
external walls are vibrating in a sawtooth way with velocity vt ■ 



dilute here and very dense in |15] . In addition, no fric- 
tion is considered in the present analysis. Very recently, 
the velocity and position fluctuations of a piston on top 
of a vibrated granular gas in presence of gravity have also 
been investigated, and a rather complex and not weW un- 
derstood behavior has been found [16 . 

We consider a cylindrical container of length Lj. divided 
into two compartments by a movable piston constrained 
to remain perpendicular to the axis of the system, taken 
as X axis. The system is sketched in Fig. [TJ The two 
compartments contain A^i and N2 smooth inelastic hard 
spheres (d = 3) or disks {d = 2), respectively. A sub- 
script i — 1,2 will be used to identify properties on each 
side of the piston. The particles within each subsystem 
are equal and their mass and diameters will be denoted 
by rrii and ai, i — 1,2, respectively. Collisions between 
particles are characterized by velocity-independent coeffi- 
cients of normal restitution a^, < < 1. On the other 
hand, collisions between particles and the piston, as well 
as collisions of particles with the walls of the container 
are taken as elastic. The piston moves without friction 
with the walls and it is modeled as having mass M and 
negligible width. Energy is continuously supplied to the 
system through the external walls. All them vibrate with 
a sawtooth velocity profile and an amplitude much smaller 
than the mean free path of the particles in their vicinity. 
That means that to describe collisions, the walls can be 
considered as fixed and with a constant velocity Vb per- 
pendicular to them, and addressed towards the interior of 
the system pj] . 

The focus here will be on the steady state position xp 
eventually reached by the piston as a function of the pa- 
rameters defining the system. This is equivalent to study 
the way in which the fixed total volume is divided into the 
two subsystems by the movable piston. In the following 
analysis, the rather drastic approximation that the granu- 
lar gas inside each of the compartments can be treated as 
homogeneous will be made. The expectation is that this 
simplification does not affect the qualitative behavior of 
the system, although it certainly will modify the accuracy 
of the quantitative predictions. Moreover, attention will 
be restricted to those situations in which both gases are 



very dilute, so that their equations of state have the form 
Pi — UiTi, where Pi is the pressure, rii the number density, 
and Tj the granular temperature of the gas in compart- 
ment i. As usual, the latter is defined from the second 
moment of the velocity distribution with the Boltzmann 
constant formally set equal to unity. It is 
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where Sp is the area (d = 3) or length {d — 2) of the 
piston and, therefore, also the section of the container. 
In the steady state, mechanical equilibrium of the piston 
requires that pi = P2 = P, i.e. 



(2) 



Although the stationarity of xp does not imply by itself 
the same for the temperatures or the pressure, only sta- 
tionary steady sates of the whole system will be addressed 
in the following. Then energy balance for the gases in each 
of the two compartments reads 



-NiQTi — QpiSp — QwiSwi = 0. 



(3) 



The first term of the left hand side of this equation is the 
Haff law, describing the homogeneous cooling of the gas 
as a consequence of the inelasticity of collisions [T^ . It in- 
volves the cooling rate Q for which an accurate expression 
is [UlIlO] 
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and rjQi{T) is the shear viscosity of a dilute molecular gas 
of hard spheres or disks. 
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The second term in Eq. ([3]) is the power going into com- 
partment i through the piston due to its velocity fiuctu- 
ations. This is the way in which the gases at both sides 
of the piston interchange energy. To describe it, it is con- 
venient to introduce a temperature parameter Tp for the 
piston defined from the second moment of its velocity fluc- 
tuations, namely 

M{V^) EE Tp, (8) 

with the angular brackets denoting statistical average. For 
nii/M <^ 1, i = 1,2, and assuming that Tp/Ti is of the 
order of unity, it is [T^[^ 
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Finally, the last term in Eq. ([3]) is the power injected into 
the subsystem i by the vibrating walls in contact with it. 
For the specific kind of vibration considered here [21], 



iWi = pvb- 
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Moreover, Swi is the total section (area or length) of the 
vibrating walls in contact with the gas in compartment 
i. This quantity is a function oi xp. The steady energy 
balance for the piston reads 

Qpi + Qp2 = 0. (11) 

Let us introduce dimensionless temperature parameters by 
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and also a = (cri+(T2)/2. In terms of the above quantities, 
Eqs. © and become 
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i = 1,2. In the above equations. 
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Equations (|T4)) and (|T5|) together with Eq. ([2]) rewritten 
as 

(20) 
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define a closed set of equations for the quantities xp, 9i, 
62, and 0p. It is worth to remind that Ci depends on 
through Swi- More specifically, it is 

Swi = Sp+-fxp, Sw2 ^ Sp +j{L^ - Xp), (21) 

with 7 being some geometrical factor with dimensions of 
(length)''"^. For instance, if the system is a circular cylin- 
der 7 — 2(7r S'p)^/^, while for a two-dimensional rectangu- 
lar system it is 7 = 2. 

Since the problem has many parameters, let us partic- 
ularize for the case in which all the particles are identical 
and there is the same number of them in each compart- 
ment, i.e. Ni = N2 = N, ai = a2 = a, Ci = fT2 = o". 



TTii = 77i2 = 777. Thcu, after some algebra, from Eqs. (jl4p . 
([T5|). and (HH), it is possible to derive a closed equation for 

(l-p)(/72 + A/7+l) =0, (22) 
2N^Lxad 
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where a — ai — a2 and = 61 = 62- Note that now it is 
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Once the density parameter p is known, the value of the 
temperature 9i can be obtained by combination of Eq. 
([T4I1 for 7 = 1 and 7 = 2 and using Eq. (fTSj) . The result is 
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Moreover, the temperature parameters are related by 

0P ^ p9i ^ {9id2)'^\ (26) 

Equation (j22p has the expected symmetric solution p — 1, 
i.e. 77i = 772, Xp = Lj./2. Also, it follows from Eqs. 
and (PO)) that 9i ^62 = 9p,so that the gases in both com- 
partments and the piston all have the same temperature 
parameter. This state exhibits the symmetry of the mate- 
rial parameters at both sides of the piston. But Eq. p2|) 
has also another solution corresponding to an asymmetric 
state, in which the steady position of the piston is not in 
the middle of the system but at a position defined by 
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The two solutions p+ and p_ — 
from the other by interchanging the two compartments. 
The existence of this asymmetric state requires that the 
parameters of the system verify the condition A > 2. It 
follows from Eq. ((20)) that the steady temperatures of the 
granular gases in the two compartments are different in 
the asymmetric state. 

To verify the accuracy of the above theoretical predic- 
tions, molecular dynamics (MD) simulations have been 
performed for a system of 2N — 200 equal inelastic hard 
disks {d = 2) enclosed in a rectangular container of size 
Lx x Ly. The latter corresponds to the section Sp used 
above. All the simulations to be reported in the follow- 
ing started from a symmetric situation with the piston at 
rest in the middle, xp{0) = Lx/2, and equal number of 
particles in both compartments. The initial velocity dis- 
tributions were Gaussian and with the same temperature. 
Finally, the values of and Ly were always chosen such 
that the initial number density was quite low. 

As predicted by the theory developed above, for small 
values of A defined in Eq. ([23)) . the piston remains os- 
cillating around the middle of the system. Nevertheless, 
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when the value of A increases above 2, the average position 
of the piston clearly moves away from the initial position 
until reaching a different steady average. In this way, the 
symmetry of the system is broken. In Fig. [21 the steady 
average position of the piston is plotted as a function of 
the dimensionless control parameter A for a system with 
Lj: = 2Ly = 133ct, so that the initial homogeneous density 



in both compartments is n 
plotted quantity is 

\2xp - 
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Different values of the coefficient of normal restitution a 
have been considered. Moreover, for each value of a, the 
values of A has been modified by changing the ratio M/ m 
between the mass of the piston and the mass of the par- 
ticles. A reasonably good agreement is found between 
theory and simulations, the former being given by the last 
equality in Eq. (I?fl) with p = p+, as given by Eq. (f^ . It 
must be taken into account that when the system is above 
but close to the bifurcation point, the position of the pis- 
ton does not oscillate around a unique steady position but 
also moves back and forth between the two possible ones, 
at both sides of the middle of the system. It is not an easy 
task to separate both motions as required to compare with 
the theory developed here, which does not consider fluc- 
tuations between the two steady positions. Similar results 
have been obtained for other values of the average num- 
ber density, as illustrated in Fig. |3]for n = 0.04cr~^. As 
expected, the quantitative accuracy of the theoretical pre- 
diction decreases as the density increases, since the model 
has been formulated in the low density limit. Let us also 
notice that for the geometry used in our simulations, when 
< Ly, \t \s X < 2 independently from the values of the 
other parameters of the system and, therefore, the theory 
predicts that there is no spontaneous symmetry breaking. 
This feature has been consistently confirmed by the sim- 
ulation results. 

It is interesting to verify at what extent the hypothesis 
on which the theory formulated above is based are verified 
and, in particular, the homogeneity of the subsystems. 
To measure the temperature and number density fields 
in the simulations, the system has been divided into cells 
parallel to the piston of width Lx/2Q. The fields have been 
measured once the system has reached the steady state, 
and the reported profiles have been averaged over time and 
also over trajectories. In Fig. |4]the steady average density 
and temperature profiles, n(x) and T{x), are shown for 
a system with = 2Ly = lOOcr, and a = 0.99, and 
M = 150m. This gives A ~ 2.1 and the prediction for 
the steady position of the piston is xp ~ 0.65Lx- It is 
observed that outside from a boundary layer at both sides 
of the piston, the hydrodynamic fields can be considered as 
uniform inside each compartment, although having quite 
different values in each one. The boundary layers are, at 
least partially, due to the oscillations of the piston. 
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Fig. 2: Bifurcation diagram showing the asymmetry of the 
steady position of the piston as a function of the dimension- 
less control parameter A defined in Eq. (|23p . The symbols are 
from molecular dynamics simulations and the solid line is the 
theoretical prediction given by Eq. ((2})) with p — p+. The 
values of the parameters are La: ~ "iLy = 133cr and 2N = 200, 
i.e. na^ = 0.02. Different values of a have been used as in- 
dicated in the insert. In addition, also different values of the 
ratio M/m have been employed to change the value of A for 
each value of a. 



Another quantitative prediction of the model are the 
values of the temperature ratios T2/T1 = ^2/^1 and 
Tp/Ti — 6p/0i. The theoretical predictions are given 
by Eqs. (PT|) and (|?7)) . They are compared with the simu- 
lation results for a system with = Ly = lOOcr in Fig. [Sj 
When the system is in an asymmetric state, the temper- 
ature of each of the compartments has been measured by 
omitting the boundary layers next to the piston. Again, 
given the approximations involved in the formulation of 
the model, the accuracy can be considered as satisfactory. 

In summary, we have found a spontaneous symmetry 
breaking in two vibrated granular gases separated by an 
adiabatic piston. Moreover, the energy balance equations 
^ and pT|) . based on simple kinetic theory arguments, 
capture surprisingly well the qualitative and quantitative 
features observed in the molecular dynamics simulations. 

Of course, a relevant question is to know which is, if 
any, the macroscopic criterion leading to equilibration of 
the granular gases in both compartments. In this sense, 
a first conclusion of the analysis here is that there is a 
macroscopic description of the steady states of the system 
in terms of the densities of both compartments and one 
temperature parameter. From them, the other two tem- 
perature parameters and the pressure can be determined, 
if needed. Also, it is possible to compare the energy dis- 
sipation in the symmetric and asymmetric steady states. 
The power Pd dissipated in collisions in a steady state is 
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3: Tiie same as in Fig. (2] but for a system with Lx = 
= lOOcr, so that no-^ = 0.04. 
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This function has a minimum for p = 1, i.e. the sym- 
metric state. This shows that the system chooses between 
the allowed states the one having a maximum energy dis- 
sipation in the granular gases. This information could be 
relevant to fromulate some extremal principle governing 
the macroscopic evolution of granular fluids. 



This research was partially supported by the Ministe- 
rio de Educacion y Ciencia (Spain) through Grant No. 
FIS2008-01339 (partially financed by FEDER funds). 



REFERENCES 



[1] 



B., Thermodynamics (Wiley and Sons, New 



Callen H. 
York) 1963. 

LiEB E. H., Physica A, 263 (1999) 491. 
Kestemont E., Van den Broeck C, and Mansour M., 
EPL, 49 (2000) 143. 

Gruber Ch., Pache S., and Lesne A., J. Stat. Phys., 
112 (2003) 1177. 

Reimann p., Phys. Rep., 361 (2002) 57. 
Gruber Ch. and Lesne A., Encyclopedia of Mathematical 
Physics, edited by Francoise J. P., Naber G., and Tsun 
T.S. (Elsevier, Amsterdam) 2006, p. 160. 
Campbell, C.S., Annu. Rev. Fluid Mech., 22 (1990) 57. 
Goldhirsch I., Annu. Rev. Fluid Mech., 35 (2003) 267. 
Brilliantov N. and Poschel T., Kinetic Theory of 
Granular Gases (Oxford University Press, Oxford) 2004. 
[10] Granular Gases: Beyond the Dilute Limit, edited by Brey 
J.J., Goldhirsch I., and Poschel T., Vol. 179 (EPJ 
ST) 2009. 

[11] Brito R., Renne M.J., and Van den Broeck c, EPL, 
70 (2005) 29. 

[12] Brey J.J. and Khalil N., Phys. Rev. E, 82 (2010) 
051301. 



[2] 
[3] 

[4] 

[5] 
[6] 



[7] 
[8] 
[9] 



1 - 



3000 



2000 



1000 



0.2 



0.4 0.1 

x/L, 



0.8 1 



Fig. 4: Dimensionless steady number density n{x)/n, n — 
N/L^Ly, (solid line and left axis) and temperature 9{x) = 
T{x)/mvl (dashed line and left axis) profiles along the x- 
direction for a system with Lx = 2Ly = 100a, a = 0.99, 
and M = 150m. The vertical line indicates the theoretical 
prediction for the average position of the piston. 



[13 
[14 

[is: 

[16 
[17 

[is; 

[19 
[20 
[21 



CosTANTiNi G., Marini bettolo Marconi U., and 
PuGLisi A., EPL, 2008 (50008) . 

PiASECKi J., Talbot J., and Viot P., J. Stat. Mech.: 
Theory Exp., 2010 (P05004) . 

Lechenault F. and Daniels K.E., Soft Matter, 2010 
(3074) . 

Brey J.J. and Ruiz-Montero M.J., Phys. Rev. E, 81 
(2010) 021304. 

McNamara S. and Luding S., Phys. Rev. E, 58 (1998) 
813. 

Haff P.K., J. Fluid Mech., 134 (1983) 401. 
GOLDSHTEIN A. AND Shapiro M., J. Fluid Mech., 282 
(1995) 75. 

VAN NoiJE T.P.C. AND Ernst M.H., Granular Matter, 
1 (1998) 57. 

Brey J.J. and Ruiz-Montero M.J., Phys. Rev. E, 79 
(2009) 031305. 



p-5 



F. Author et al. 



T2/T1 
Tp/Ti 



+ 



2 



Fig. 5: Ratio of the granular temperatures of the gases in both 
compartments T2/T1 (solid line and plus symbols) and between 
the temperature parameter of the piston and that of compart- 
ment 1, Tp/Ti (dashed line and crosses) as a function of the 
control parameter A, defined in Eq. (|23[) . The lines are the the- 
oretical predictions derived here while the symbols are molec- 
ular dynamics simulation results. The size of the system is 
La. = 2L„ = 100a. 
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